Prom Linear to Nonlinear Response in Spin Glasses: 
Importance of Mean-Field-Theory Predictions 
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Deviations from spin-glass linear response in a single crystal Cu:Mn 1.5 at % are studied for a wide 
range of changes in magnetic field, AH. Three quantities, the difference TRM — (MFC — ZFC), 
the effective waiting time, tiff, and the difference TRM(t w ) — TRM(t w = 0) are examined in 
our analysis. Three regimes of spin-glass behavior are observed as AH increases. Lines in the 
(T, AH) plane, corresponding to "weak" and "strong" violations of linear response under a change 
in magnetic field, are shown to have the same functional form as the de Almeida- Thouless critical 
line. Our results demonstrate the existence of a fundamental link between static and dynamic 
properties of spin glasses, predicted by the mean-field theory of aging phenomena. 

PACS numbers: 75.50.Lk,75.40.Gb 



I. INTRODUCTION 

The properties of spin glasses in a magnetic field have 
been the subject of considerable attention. It is widely 
recognized that the Parisi replica-symmetry breaking 
ansatzEJ provides an essentially correct mean-field so- 
lution |-|for the infinite-range Shcrrington-Kirkpatrick 
model.H The spin-glass phase in this model is separated 
from the paramagnetic phase in the (T, H ) plane by the 
de Almeida- Thouless (AT) critical line.13 The situation 
is less clear in the case of finite-dimensional short-range 
models. Rigorous theoretical results show viability of t 
mean-field approach for the description of these systems,! 
but progress is impeded by significant analytical difficul- 
ties. Numerical studies have repeatedly suggested the 
existence otajji AT-type critical line at d = 3 and higher 
dimensions. B'tl However, magnetic field effects present a 
difficult challenge for computer simulations. Finite size, 
and difficulties in equilibrating large samples, do not yet 
allow a clear distinction .between the mean-field picture 
and the droplet scenario .13 Because of this, the existence 
of the spin-glass phase transition in a magnetic field is 
considered a most relevant open problem.H 

The experimental evidence for an AT-type critical be- 
havior appears ambiguous. Many real spin glasses have 
mean- field-like phase diagrams Jj with an onset of strong 
MFC — ZFC irreversibility along a certain line. This 
irreversibility is usually interpreted as a sign of replica- 
symmetry breaking. However, real spin glasses are al- 
ways out of equilibrium, and the measured AT lines are 
time dependent. Therefore, it appears impossible to ob- 
tain information about the equilibrium phase diagram di- 
rectly from experiments. We show in this paper, however, 
that the study of magnetic field effects on the nonequi- 
librium dynamics of spin glasses can shed light on the 
magnetic properties of the equilibrium spin-glass state. 

Our motivation for this work is twofold. First, effects of 
magnetic field are understood fairly well within the mean- 
field theory, where they are derived from first principles. 
The minimum possible overlap, q m in{H), of two states in 
the presence of a magnetic magnetic field, if, plays an 



important role in this approach. Yet, as usually happens 
in spin-glass physics, the theoretical predictions cannot 
be easily related to experimentally observable phenom- 
ena. Even those models of spin-glass dynamics that are 
based on the mean-field-like hierarchical picture of phase 
space tend to treat magnetic field effects phenomenolog- 
ically in terms of the Zeeman energy. In this paper, we 
show that experimental results on the magnetic field de- 
pendence of spin-glass dynamics are consistent with the 
mean-field description. 

The second reason for this work is practical. Spin-glass 
relaxation properties are usually studied under a change 
in magnetic field. It is often believed that the subsequent 
response is linear in magnetic field "for reasonably small 
field values (say < 10 G)" .□ We show in this paper that 
the very definition of what one means by "reasonably 
small" fields is impossible without knowledge of the spin- 
glass phase diagram. This knowledge becomes vital when 
results for different temperatures or different samples are 
compared. 

The paper is organized as follows. The next Sec- 
tion describes a theoretical picture underlying our anal- 
ysis. Sec. III. A presents results of magnetization mea- 
surements. Sec. III.B is devoted to experimental results 
for the effective waiting time. In Sec. III.C, the waiting- 
time dependence of the measured quantities is discussed 
in detail. Sec. IV summarizes our conclusions. 



II. THEORETICAL BACKGROUND 

A. Mean-field dynamics 

The main obstacle to an experimental test of mean- 
field-theory predictions is the problem of relating ex- 
perimentally accessible spin-glass dynamics to the static 
equilibrium properties of the spin-glass state, described 
by the Parisi solution. Phenomenological phase-space 
models have provided important insights intp ..possible 
physical mechanisms for spin-glass relaxation.tLHltil A re- 
cent theoretical breakthrough has been achieved within 
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eapz_.of paging phenomena, proposed by 
.ESOaa Because these theoretical de- 



the mean-field the 
Cugliandolo et all 
velopments have profound significance for the interpreta- 
tion of our experimental results, we shall briefly review 
them here. 

The j— - well-known fluctuation-dissipation theorem 
(FDT)IIa establishes a link between the linear response 
of a system to an external perturbation, and the 
fluctuation properties of the system in thermal equi- 
librium. The response of a spin system at time t to 
an instantaneous field at time t', R(t,t'), and the spin 
autocorrelation function, C(t,t'), are defined as follows: 



N 



R(t,t') = (l/N)Y,HS i (t))/5h(t') 



i=l 



N 



C(t,t') = (l/N)^2{S i (t)S i (t')) 



(1) 



(2) 



In thermodynamic equilibrium, both functions are time- 
translation invariant, and related by the fluctuation- 
dissipation theorem ((3 = 1/fcsT): 



R eq (t-t')=pdC eq (t-t')/dt' 



(3) 



The violation of this theorem in a general off-cquih^rium 
situation is described by the function X(t, t') < 11BBK3 



R(t,t') = /3X(t,t')dC(t,t')/dt' 



(4) 



Violation of the fluctuation-dissipation theorem is asso- 
ciated with a deviation from linear response. A second- 
order nonlinearity does not appear because the magne- 
tization changes sign when the field is reversed. It has 
been showrcZl that the presence of a third-order nonlinear 
response turns the fluctuation-dissipation theorem into a 
fluctuation-dissipation inequality. This is consistent with 
the fact that X(t,t') in Eq. (4) is generally less than 
unity. 

When these formulas are applied to spin-glass dynam- 
ics, the relevant times are t = t w + r and t' = t w , where 
t w and r are the waiting time and observation time after 
the waiting time, respectively. The dynamical quantities 
C and X can be related to their equilibrium counter- 
parts if both t w and r are sent to infinity. However, 
time-translational invariance does not hold off equilib- 
rium, and the limiting values depend on how the luni|te. 
are taken. One can write for the correlation functionilS'Ell 



lim lim C(t w + t, t w ) = q E A , (5) 

r — >oc t w — >oo 



and 



lim C(t w +r,t w ) = q n 



(6) 



Eq. (5) is a dynamical definition of the Edwards- 
Anderson order parameter, qea, as the value of the corre- 
lation function at the limit of validity of the fluctuation- 
dissipation theorem. Eq. (6) expresses the property of 



weak ergodicity breaking in spin glasses. If the wait- 
ing time is finite and the magnetic field is zero (giving 
q m in — 0), the system is able to escape arbitr.arily far 
from the configuration it had reached at t — t w x3 

The fluctuation-dissipation ratio X(q) is defined as the 
limit of the function X(t,t'), taken along a path in the 
{t,t') plane, characterized-by a given value q of the cor- 
relation function C(t,t')£3 The most important result 
of the theory is the_fact that, under the assumption of 
stochastic stability ,li3 X(q) is equal to the equilibrium 
order parameter function x(q): 



X(q) = lim 



R(t,t') 



t,t'^oo; c= g (3dC{t,t')/dt' 



(7) 



The function ff(q) is an integral of the Parisi order pa- 
rameter P(q)p which can be introduced for short-range 
spin glasses within the standard replica-symmetry break- 
ing formalism^ According to Eq. (7), there is a deep 
relationship between spin-glass dynamics and the static 
equilibrium properties of the spin-glass state. It implies 
that "in any finite-dimensional system, replica-symmetry 
breaking and aging in the response functions either ap- 
pear together or do not appear at all" .E3 

We now consider the physical (experimental) signifi- 
cance of these theoretical results. It is well known that 
there are two main regimes of spin-glass relaxation. In 
the equilibrium, or stationary regime, observation times 
r < t,,, the fluctuation-dissipation theorem holds and 
the response function depends only on r. In the nonequi- 
librium, or aging regime, observation times r > t w , the 
FDT is violated and the relaxation depends on t w for au$4 
t. Both regimes have bcen-observed experimentallyll3 
and studied numerically.liijEj The transition from one 
regime to the other is marked by a peak in the relax- 
ation rate, corresponding to strong violation of FDT at 

T ~ t w . 

These effects can be naturally explained within the 
mean-ficld—thcory of aging phenomena. It has been 
suggcstccH that the function X(t,t'), introduced to de- 
scribe violations of the FDT in Eq. (4), depends on its 
time arguments only through the correlation function: 
X = X[C(t, t')]. In the equilibrium regime, the correla- 
tion function decreases rapidly (on the linear time scale) 
from 1 to qEA, and X(C) = 1. In the aging regime, 
the correlation, fun ction relaxes slowly from qEA to q m ini 
and X{C) < l.EillH3 The transition from one regime to the 
other is a direct consequence of replica-symmetry break- 
ing at q < qEA- 



B. Chaotic nature of the spin-glass state 

Another important issue is the chaotic nature of the 
spin glass state with respect to magnetic field. It has 
been demonstrated numerically that a small change in 
external field leads to a considerable reorganization of 
a spin configuration. E3 The Parisi solution suggests that 
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an average equilibrium overlap between two states at dif- 
ferent but similar magnetic fields, (hi — h 2 ) 2 N j$4 ljj is 
equal to gjnin j i- e - the minimum possible overlapBil Anal- 
ysis of fluctuations around the Parisi solution, carried out 
by Kondor,E3 demonstrates that the correlation overlap 
function Cjj(r) for two spins, i and j, at a distance r, 
behaves as 



C H (r) = (S l S J ) H (S l S J ) cx exp(-r/&r) . (8) 

This means that the projection of the correlation (SiSj)H 
at field H onto the correlation (SiSj)o at zero field van- 
ishes beyond the finite characteristic length £#. Near 
T g and at low fields, this magnetic correlation length 
is simply = \/q min (H). In the mean-field theory, 
q m in(H) oc H 2 ^, so that £h diverges rapidly as the field 
goes to zero. 

The behavior of the correlation function, Eq. (8), is 
a consequence of replica-symmetry breaking. The low- 
temperature spin-glass phase has an essentially infinite 
number of pure equilibrium states. Each of them is char- 
acterized by an infinite correlation length. These states 
have equal free energies per site, except for differences of 
the order 0(1 /N). Only a few states with the lowest en- 
ergies contribute significantly to the partition function. 
Because of the small energy differences, any small (but 
finite) amount of energy, added to the system, is enough 
to reshuffle the Boltzmann weights of the different states 
and thus completely reorganize the equilibrium spin con- 
figuration. Application of a magnetic field is an example 
of such a perturbation. 

The magnetic correlation length, has the following 
meaningful When a magnetic field H is applied to the 
system, the minimum possible overlap of two states is 
equal to q m i n (H). Consequently, all the states having 
overlaps q < q m in are suppressed by the field. Their free 
energies increase, and they acquire the finite correlation 
length The spatial spin correlations, corresponding 
to these states, survive only within this range. All the 
other pure states are still characterized by an infinite 
correlation length. At the AT line, where q m in — qEA, all 
the states have a finite correlation length, and the system 
becomes paramagnetic. Thus, a change in magnetic field 
has a randomizing effect on the spin-glass state. The 
ratio qmin/ qEA is a natural measure of this effect. 

The chaotic nature of the spin-glass state is also re- 
flected in the phenomenological droplet model. The 
magnetic correlation length is determined as an average 
droplet size for which the Zeeman energy is equal to the 
energy of the droplet excitation. It is given by the fol- 
lowing expression£3 

i H oc tf- 2 /( d - 29 ) . (9) 

Even though this result, with d = 3 and 8 « 0.2, is sim- 
ilar to the result of the mean-field theory, the physics is 
quite different. Spin-glass properties in this model are 
governed by low-energy excitations of the ground state, 
created by coherent flipping of compact clusters of spins. 



It is suggested that the magnetic field, H, would flip 
all the droplets with sizes greater than and thus de- 
stroy all spin-glass correlations beyond this length scale. 
Therefore, at any nonzero field, there is only a paramag- 
netic state with the finite correlation length 

In both these approaches, the chaotic nature of the 
spin-glass state with respect to magnetic field leads to 
departures from linear response as AH increases. 



III. EXPERIMENTAL RESULTS AND 
ANALYSIS 

The purpose of this paper is the study of violations of 
the fluctuation-dissipation theorem as a function of mag- 
netic field change AH. According to Eq. (7), aging dy- 
namics of the spin-glass state is ultimately determined by 
static equilibrium properties. Any violation of the FDT 
contains information about replica-symmetry breaking. 
Therefore, studies of the gradual deviation from linear 
response as AH increases can provide insight into the 
nature of the spin-glass phase diagram. Direct experi- 
mental determination of the correlation function, Eq. (2), 
requires sophisticated measurements of the time depen- 
dent magnetic noise spectrum. E3 Instead, we make use of 
quantities which can be obtained from magnetic suscepti- 
bility measurements. The first quantity of interest is the 
difference of the values of remanence, measured in TRM 
and ZFC experiments: TRM - (MFC - ZFC). It has 
been shownE3 that, within the linear regime, this quan- 
tity is zero, provided that all three magnetizations are 
measured at the same time after the initial quench. The 
second quantity of interest is the effective waiting time, 
tff , defined as the value of the observation time r where 
the relaxation rate S(r) = -dTRMfcr, t w , AH)/d\ogT 
has a maximum. It has been arguedEJ that, within the 
linear response regime, tiff w t w , so that the peak in 
S(t) is essentially unaffected by a small magnetic field 
change. The third quantity we study experimentally is 
the difference TRM(t w ) - TRM (t w = 0). It describes a 
change in the measured magnetization as a result of the 
waiting time, and thus allows closer examination of aging 
phenomena. 

All experiments were performed on a single crystal of 
Cu:Mn 1.5 at %, a typical Heisenberg spin glass with a 
glass temperature of about 15.2 K. Results of conven- 
tional studies of the phase diagram for this sample will 
be reported elsewhere. A commercial Quantum Design 
SQUID magnetometer was used for all the measurements. 



A. Measurements of TRM, MFC, and ZFC 

Fig. 1 exhibits the dependence oiTRM, MFC-ZFC, 
and their difference, on the field change AH , measured 
at T = 12.0 K. All data points were taken at the same 
short observation time, r f=a 40 s, with zero waiting time, 
and an effective cooling time of about 600 s. Error bars 
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AH (kOe) T (K) 



FIG. 1: The values of the remanent magnetic moment, mea- 
sured at T = 12.0 K = 0.79 T g and t w = 0. The magnetic field 
is changed from H = AH to H = (TRM experiment) and 
from to AH (ZFC experiment). The two arrows indicate the 
critical field changes, A#i w 120 Oe and AH 2 « 420 Oe, for 
the weak and strong linear response violations, respectively. 



are smaller than the symbol sizes. The same will apply to 
all figures without error bars. Fig. 1 shows that three dif- 
ferent types of spin-glass behavior can be distinguished 
for different ranges of magnetic field change, AH. For 
field changes from AH = to AH X « 120 Oe, the dif- 
ference TRM - (MFC - ZFC) is zero to within our ex- 
perimental accuracy. Between AH\ and AH2 « 420 Oe, 
there are weak deviations from linear response. At AH = 
AH2, corresponding approximately to the peak in the 
MFC — ZFC remanence, violation of linear response be- 
comes strong, and the difference TRM - (MFC - ZFC) 
is a linear function of AH with a large slope. 

The temperature dependences of the critical field 
changes AHi and AH2 are presented in Fig. 2 and Fig. 3. 
The figures also show the critical AT line, determined 
for the same sample at the same observation and cooling 
times, using the onset of strong MFC —ZFC irreversibil- 
ity as the signature of the spin-glass phase transition. We 
assume that this dynamical line approximates the equi- 
librium AT line for infinite waiting time. The difference 
MFC - ZFC in Fig. 1 is not zero at the AT field of 
w 1400 Oe because transverse freezing above the AT line 
leads to a weak longitudinal irreversibility.tB 

One sees from Fig. 2 that all three lines in the (T, AH) 
plane have the same functional form, T g — T. = a AH 2 / 3 , 
typical of the equilibrium AT critical linetj in the SK 
model. Moreover, the two linear-response violation lines, 
AHi(T) and AH2(T), determine dynamical transition 
temperatures which are very close to the actual glass 
temperature T g . It would thus be incorrect to say that 
the critical AT-type line plays a role at high fields only. 
Our results demonstrate that this line manifests itself dy- 
namically even at very low fields. A critical field change 



FIG. 2: The dynamical AT line AH A t(T), the strong linear 
response violation line AH 2 (T), and the weak linear response 
violation line AHi(T), determined from the linear fits to ex- 
perimental data, plotted as AH 2 ^ 3 vs. T. 
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T (K) 

FIG. 3: The same lines as in Fig. 2, but plotted as AH(T). 
Note that the two linear response violation lines below the AT 
line represent a dynamic crossover from one relaxation regime 
to the other. They are not features of an equilibrium phase 
diagram. 

for a given degree of linear response violation seems to 
be a constant fraction of the AT field. It means, from a 
practical point of view, that results obtained at different 
temperatures can be directly compared only if they have 
the same AH / Hat{T) ratio. 

A comment should be made at this point. It has been 
arguedo that a T cx H 2 ! 3 dependence for transition lines 
is not a unique feature of the mean-field theory. A similar 
power law could be obtained from purely dynamical con- 
siderations without using the concept of the spin-glass 
phase transition^ A dynamical freezing line, Tf(H,u), 
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FIG. 4: The function TRM - {MFC - ZFC), describing 
violation of linear response, 3,t tyj — 30 min and t w = 0. The 
second and the third arrows from the left indicate, respec- 
tively, the larger value of AHi and the lower value of AH2 
for the longer waiting time. 

appears in the droplet model as well.0 In the present 
analysis, we do not rely on the existence of the H 2 I Z de- 
pendence by itself. Our argument is based on the fact 
that the experimental AT and linear response violation 
lines have the same functional form, as illustrated in 
Fig. 2 and Fig. 3. This result suggests that there is a 
close relationship between the static and dynamic prop- 
erties of the spin-glass state, assuming that the measured 
AT line is related to the equilibrium one. 

Of course, measured values of the critical field changes 
depend on both the waiting time and the observation 
time. This point is illustrated in Fig. 4 where values 
of TRM - {MFC - ZFC) are plotted for t w = and 
t w = 30 min. As the waiting time becomes larger, AHi 
increases, while AH2 decreases. Thus, longer equilibra- 
tion times make the system less susceptible to external 
perturbations, but only up to a certain point. This issue 
will be discussed further in Sec. III.C. 



B. Measurements of the effective waiting time 

Let us now turn to a discussion of the effective waiting 
time, t^/ . Fig. 5 exhibits the dependence of log 10 (tj« ) 
on the magnetic field change AH for three different wait- 
ing times: t w = 30 min, 5 min, and 0. The temperature 
is the same as in Fig. 1. One can easily discern three 
regimes of FDT violation, corresponding to the three 
regimes of deviation from linear response seen in Fig. 1. 

As the magnetic field change, AH, increases from zero 
up to about 100 Oe, the effective waiting time decreases 
only slightly. This means that FDT holds for observation 
times less than tffi ~ t w . Of course, if the waiting time 
is short, tffi is determined by the cooling procedure. 



FIG. 5: Logarithm of the effective waiting time, t%, as a 
function of AH, measured at T = 12.0 K = 0.79 Tg for three 
waiting times. The dotted lines are linear fits in the interval 
100. ..250 Oe. They are only guides for the eye. 



In the interval from 100 Oe to about 400 Oe, the effec- 
tive waiting time drops sharply, and the field dependence 
of its logarithm is linear. In this case, nonequilibrium be- 
havior appears much earlier, and the FDT holds only at 
short times r <C tiffi <C t w . This corresponds to the weak 
deviation from linear response between AH± and AH2 in 
Fig. 1. 

When the field change exceeds 400 Oe, the FDT is 
strongly violated. One can see from Fig. 5 that, in this 
regime, the curves for long waiting times approach the 
curve for zero waiting time, and the waiting time depen- 
dence almost disappears. This situation corresponds to 
the strong deviation from linear response at short obser- 
vation times above AH2 in Fig. 1. Therefore, the break 
from the linear dependence of log 10 (i^) on AH, ob- 
served at large field changes, is directly related to strong 
nonlinearity in the spin-glass response. The effective 
waiting time in this case is determined primarily by the 
zero waiting time results, dependent on the experimental 
protocol. 

The waiting time dependence in Fig. 5 supports the 
conclusion, drawn from Fig. 4, that an increase in waiting 
time leads to an expansion of the linear response region to 
higher values of AH. An increase in the observation time 
t has the opposite effect. The field dependence is more 
pronounced at longer times. This accounts for the fact 
that the critical field changes, extracted from the effec- 
tive waiting time experiments (Fig. 5), seem to be lower 
than those obtained from the short-time magnetization 
measurements (Fig. 4). Another interesting feature of 
our results is that the dependence of log 10 (f^ ') on AH 
for t w = 0, determined exclusively by the cooling proce- 
dure, does not exhibit the low-field plateau seen at longer 
waiting times. This behavior is analyzed in Sec. III.C. 

The temperature dependence of the effective waiting 
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Logarithm of the effective waiting time, t%, mea- 



sured as a function of AH for t v 



30 min at four tem- 



peratures: Tr = T/T g = 0.70, 0.75, 0.79 and 0.85. Note the 
similarity in shape and pronounced difference in the magnetic 
field scales. 



FIG. 8: Comparison of TRM(AH) curves measured at T — 
12.0 K — 0.79 T g for different waiting times. The curve for 
t w — is referred to as ZTRM. As t w increases, the peak 
in TRM shifts to the left, while the peak in TRM - ZTRM 
moves to the right. 




2.0-| , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

AH /AH 2 

FIG. 7: The experimental results of Fig. 6, plotted vs. 
AH/AH2. The critical field change AH2, corresponding to 
strong violation of linear response, is proportional to the AT 
field, Hat{T). The straight line is a guide for the eye. 



time is exhibited in Fig. 6 and Fig. 7, where log 10 (i^) 
vs. AH data for t w = 30 min are plotted for four dif- 
ferent temperatures. The results for the smallest field 
change of 10 Oe are not exactly the same because the 
effective cooling time increases as the measurement tem- 
perature is lowered. It is quite evident from Fig. 7, how- 
ever, that the four curves scale rather well together if 
plotted vs. AH/AH 2 . This means that the critical AT 
line sets a characteristic magnetic field scale at any tem- 
perature below T g , and that it plays an important role 
for all aspects of spin-glass dynamics. 



C. Waiting time dependence of AH\ and AH2 

The waiting time dependences of the characteristic 
field changes, corresponding to the weak and strong lin- 
ear response violations, deserve special attention. Our 
experimental results, exhibited in Fig. 4 and Fig. 5, sug- 
gest that AHi increases with t w , while AH 2 diminishes. 

In order to check this conclusion, we have measured 
the field dependence of the thermorcmanent magnetiza- 
tion, TRM(r,t w , AH), for three different waiting times: 
t w = 30 min, 5 min, and 0. The measurement temper- 
ature is T = 12.0 K, and the observation time after the 
field change is r = 40 s. For convenience, we shall refer 
to the zero- waiting-time dependence as ZTRM(t, AH). 
This nonequilibrium relaxation function is determined by 
the cooling process and not by the waiting time. The 
results are exhibited in Fig. 8. One can sec that the 
maxima in the TRM curves shift towards lower fields as 
the waiting time increases. To examine evolution of the 
measured magnetization curves due to aging phenomena, 
we consider the differences TRM — ZTRM, where both 
TRM and ZTRM are measured at the same observation 
time. Fig. 8 suggests that the maxima in TRM — ZTRM 
correspond to higher fields at longer waiting times. 

Fig. 9 exhibits the logarithm of the ratio {TRM — 
ZTRM) I AH. Two different types of behavior can be 
clearly distinguished. The ratio increases as long as the 
field change is less than about 120 Oe, and decreases 
at higher AH. The field dependence of its logarithm 
can be approximated by straight lines in both regimes. 
The results in Fig. 9 can be understood if we compare 
them with the data in Fig. 5. At low field changes, 
the effective waiting time for ZTRM decreases more 
steeply with increasing AH than the effective waiting 
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FIG. 9: Logarithm of the normalized increase in the ther- 
moremanent magnetization, (TRM— ZTRM) /AH, as a func- 
tion of AH at a constant observation time r = 40 s. The 
dashed lines are linear fits in the interval 0...100 Oe. The 
solid lines are linear fits in the interval 200. ..1000 Oe. The 
arrows indicate points where nonlinearity appears. 



time for TRM, which has a plateau in this region. The 
decay of the ZTRM becomes faster in time, and the 
ratio (TRM - ZTRM)/ AH, measured at fixed obser- 
vation time, increases with AH. At higher magnetic 
field changes the drop in the effective waiting time for 
the TRM is greater than for the ZTRM. As a result, 
the ratio (TRM - ZTRM) /AH decreases. According 
to Fig. 9, the transition from one regime to the other oc- 
curs at a higher magnetic field change for a longer waiting 
time. This observation supports our conclusion that the 
plateau in the field dependence of the effective waiting 
time broadens as the system equilibrates. 

Fig. 10 provides further insight into the nature of this 
effect. It displays the same quantity, the logarithm of 
(TRM - ZTRM)/ AH, as a function of time at constant 
AH = 10 Oe. The experimental protocol is the same 
as before, except for the different waiting times. One 
can see a deep similarity between the results in Fig. 9 
and Fig. 10. The waiting time dependence exhibited in 
Fig. 10 is readily understandable: the system stays longer 
in the quasiequilibrium regime for longer t w , so the max- 
imum in TRM — ZTRM shifts towards longer observa- 
tion times. Comparison with Fig. 9 suggests that the 
same is true for the field dependence: the system in the 
quasiequilibrium regime can sustain a stronger perturba- 
tion for a longer waiting time. 

This similarity between the effects of the observation 
time and of the field change on the spin-glass state 
is a natural consequence of replica-symmetry breaking. 
The weak ergodicity breaking scenario, mentioned in 
Sec. II. A, suggests that, as the observation time increases, 
the system can evolve very far from the state at time 
t w , with the minimum correlation set by q m in- This re- 



FIG. 10: Logarithm of the normalized increase in the ther- 
moremanent magnetization, (TRM —ZTRM) / AH , as a func- 
tion of log(r) at a constant field change AH — 10 Oe. The 
straight lines are guides for the eye. The arrows indicate tran- 
sitional points. 



suits from the fact that the average relaxation time is 
infinite because of the exponential distribution of free 
energies of different states. Chaotic nature of the spin- 
glass state with respect to magnetic field, discussed in 
Sec. II. B, forces the new equilibrium configuration after 
a field change to have the minimum overlap q m %n with 
the old configuration. This phenomenon is also related 
to existence of an essentially infinite number of states 
with very similar free energies. Thus, both an increasing 
observation time and an increasing field change make the 
spin-glass state less correlated with the state at time t w . 

It is of interest to examine the waiting time depen- 
dences of AH\ and AH2 for the two "competing" de- 
scriptions of spin-glass dynamics: the droplet model and 
the phase-space picture^ 

In the droplet model,EJ if R(t w ) is an average droplet 
size, and L T is an observation length scale, two limiting 
cases can be distinguished. For R(t w ) <C the FDT is 
violated when L T ~ R(t w ), and the field does not play a 
significant role. For £ H <C R(t w ), the FDT is violated as 
soon as L T ~ and the waiting time is relatively unim- 
portant. These regimes correspond approximately to ex- 
perimental regimes with AH < AH\ and AH > AH2, 
as discussed above. There is also an intermediate regime 
with R(t w ) ~ £h. This regime is very interesting phys- 
ically, because violation of the FDT depends on inter- 
play among all three length scales: L T ~ R(t w ) ~ 
Unfortunately, no predictions for this regime are given 
within the droplet modelB3 It has been arguedEI that, if 
the crossover between linear and nonlinear regimes is de- 
fined by a condition R(t w ) = the field change, needed 
to provoke nonlinear relaxation, should decrease with in- 
creasing t w . This argument can explain the waiting time 
dependence of AH2- However, it fails in the case of AH\, 
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which increases with t w . Therefore, the experimentally 
observed waiting time dependence of AHi appears rather 
counterintuitive within the droplet scenario. 

The phase-space picture of spin-glass dynamics does 
provide a consistent explanation for these phenomena. In 
this approach, evolution of a system in the phase space 
can be viewed as a series of transitions among traps, 
separated by free-energy barriers.LU As the system equi- 
librates, it encounters traps with higher barriers, and 
these traps increasingly xesemble the pure states, con- 
tributing to equilibrium £3 The quasiequilibrium relax- 
ation regime at short observation times corresponds to 
evolution within each trap, and the subsequent nonequi- 
librium behavior is related to evolution from trap to trap. 
As the waiting time, t w , increases, the system has to over- 
come a higher effective barrier to leave a trap. This takes 
a longer time, or (if the observation time is fixed) a larger 
field change. Therefore, the characteristic field change 
for the weak linear response violation, AH\, increases 
with t w . However, after overcoming a higher barrier, the 
system can explore a broader free-energy landscape due 
to hierarchical structure of the phase space. It takes less 
time or a smaller field change to produce a state very dif- 
ferent from the one at t w . The nonequilibrium behavior 
is thus more affected by external perturbations for longer 
waiting times. Therefore, the characteristic field change 
for the strong linear response violation, AH2, decreases 
with t w . 

The transition from the quasiequilibrium to the 
nonequilibrium regime is better defined after longer 
waiting times, both as a function of the observation time 
and as a function of the field change. This conclusion 
agrees with results of numerical simulations, which 
show the same effect when susceptibility is studied as a 
function of correlation.ES 



IV. SUMMARY 



The usual method for studying the spin-glass phase 
diagram is the observation of the MFC — ZFC ir- 
reversibility as a function of temperature at a fixed 
magnetic field, H = AH. A different approach is 
employed in the present paper. The temperature is 
kept constant and " irreversibility of the irreversibility" , 
TRM - (MFC - ZFC), is measured as a function of 
magnetic field change AH. Increasing AH leads to a 
gradual deviation from linear response, and, equivalently, 
to a violation of the fluctuation-dissipation theorem. 
This violation, according to the mean-field theory of 
aging phenomena, is directly related to replica-symmetry 
breaking. Our experiments show that both methods give 
the same functional form for the critical lines, suggesting 
a mean-field-like phase diagram. This conclusion is 
further supported by measurements of the effective 
waiting time, ifff . The study of waiting time effects 
on the linear response violation suggests the validity of 
the phase-space picture for spin-glass dynamics. Our 
results demonstrate the existence of a fundamental link 
between static and dynamic properties of spin glasses, 
predicted by the mean-field theory of aging phenomena. 



We would like to thank Professor J. A. Mydosh for pro- 
viding us with the Cu:Mn single crystal sample, prepared 
in Kamerlingh Onnes Laboratory (Leiden, The Nether- 
lands). We are also grateful to Dr. G. G. Kenning for 
numerous interesting discussions and to Dr. J. Hammann 
and Dr. E. Vincent from CEA Saclay (France) for help 
and advice. 



M. Mezard, G. Parisi, and M. A. Virasoro, Spin Glass The- 
ory and Beyond (World Scientific, Singapore, 1987). 

D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 
1792 (1975). 

J. R. L. de Almeida and D. J. Thouless, J. Phys. A 11, 
983 (1978). 

E. Marinari, G. Parisi, F. Ricci-Tersenghi, J. J. Ruiz- 
Lorenzo, and F. Zuliani, J. Stat. Phys. 98, 973 (2000). 

S. Caracciolo, G. Parisi, S. Patarnello, and N. Sourlas, Eu- 
rophys. Lett. 11, 783 (1990); E. R. Grannan and R. E. Het- 
zel, Phys. Rev. Lett. 67, 907 (1991); J. C. Ciria, G. Parisi, 

F. Ritort, and J. J. Ruiz-Lorenzo, J. Phys. I (France) 3, 
2207 (1993). 

E. Marinari, G. Parisi, and F. Zuliani, J. Phys. A 31, 1181 
(1998). 

F. Krzaka la, J. Houdayer, E. Marinari, O. C. Martin, and 

G. Parisi, |cond-mat/0107366. 



G. G. Kenning, D. Chu, and R. Orbach, Phys. Rev. 
Lett. 66, 2923 (1991); F. Lefloch, J. Hammann, M. Ocio, 
and E. Vincent, Physica B 203, 63 (1994) and references 



therein. 

9 G. J. M. Koper and H. J. Hilhorst, J. Phys. (France) 49, 
429 (1988). 

10 M. Lederman, R. Orbach, J. Hammann, M. Ocio, and 
E. Vincent, Phys. Rev. B 44, 7403 (1991); J. Hammann, 
M. Lederman, M. Ocio, R. Orbach, and E. Vincent, Phys- 
ica A 185, 278 (1992); Y. G. Joh, R. Orbach, and J. Ham- 
mann, Phys. Rev. Lett. 77, 4648 (1996). 

11 J.-P. Bouchaud, J. Phys. I (France) 2, 1705 (1992); J.- 
P. Bouchaud and D. S. Dean, J. Phys. I (France) 5, 265 
(1995); E. Vincent, J.-P. Bouchaud, D. S. Dean, and 
J. Hammann, Phys. Rev. B 52, 1050 (1995). 

12 L. F. Cugliandolo and J. Kurchan, Phys. Rev. Lett. 71, 
173 (1993); J. Phys. A 27, 5749 (1994); Phil. Mag. B 71, 
501 (1995). 

13 S. Franz and H. Rieger, J. Stat. Phys. 79,749 (1995). 

14 E. Marinari, G. Parisi, F. Ricci-Tersenghi, and J. J. Ruiz- 
Lorenzo, J. Phys. A 31, 2611 (1998). 

15 S. Franz, M. Mezard, G. Parisi, and L. Peliti, Phys. Rev. 
Lett. 81, 1758 (1998); J. Stat. Phys. 97, 459 (1999). 



9 



R. Kubo, Rept. Progr. Phys. 29, 255 (1966). 

L. P. Levy and A. T. Ogielski, J. Math. Phys. 30, 683 

(1989). 

J.-P. Bouchaud, L. F. Cugliandolo, J. Kurchan, and 

M. Mezard in Spin Glasses and random fields ed. by 

A. P. Young (World Scientific, Singapore, 1997). 

M. Ocio, H. Bouchiat, and P. Monod, J. Phys. Lett. 

(France) 46, L647 (1985); L. Lundgren, P. Nordblad, 

P. Svedlindh, and O. Beckman, J. Appl. Phys. 57, 3371 

(1985). 

J. O. Andersson, J. Mattsson, and P. Svedlindh, Phys. Rev. 
B 46, 8297 (1992). 

L. F. Cugliandolo, D. R. Grempel, J. Kurchan, and E. Vin- 
cent, Europhys. Lett. 48, 699 (1999). 
G. Parisi, Nuovo Cimento D 20, 1221 (1998). 

F. T. Bantilan and R. G. Palmer, J. Phys. F 11, 261 (1981). 

G. Parisi, Phys. Rev. Lett. 50, 1946 (1983); Physica A 
124, 523 (1984). 



I. Kondor, J. Phys. A 22, L163 (1989). 

F. Ritort, Phys. Rev. B 50, 6844 (1994); Phil. Mag. B 71, 

515 (1995). 

D. S. Fisher and D. A. Huse, Phys. Rev. B 38, 373 (1988); 
38, 386 (1988). 

H. Bouchiat and M. Ocio, Comments Cond. Matt. Phys. 
14, 163 (1988). 

L. Lundgren, P. Nordblad, and L. Sandlund, Europhys. 
Lett. 1, 529 (1986). 

C. Djurberg, J. Mattsson, and P. Nordblad, Europhys. 
Lett. 29, 163 (1995). 

M. Gabay and G. Toulouse, Phys. Rev. Lett. 47, 201 
(1981); D. M. Cragg, D. Sherrington, and M. Gabay, Phys. 
Rev. Lett. 49, 158 (1982). 

L. E. Wenger and J. A. Mydosh, Phys. Rev. B 29, 4156 
(1984). 



